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ABSTRACT 


A  two-equation  turbulence  model  of  the  K  —  e  type  was  recently  derived  by  Yakhot  & 
Orszag  based  on  Renormalization  Group  (RNG)  methods.  It  was  later  reported  that  this 
RNG  based  model  yields  substantially  better  predictions  than  the  standard  K  —  e  model  for 
turbulent  flow  over  a  backward  facing  step  -  a  standard  test  case  used  to  benchmark  the 
performance  of  turbulence  models  in  separated  flows.  The  apparent  improvements  obtained 
from  the  RNG  K  —  e  model  were  attributed  to  the  better  treatment  of  near  wall  turbulence 
effects.  In  contrast  to  these  earlier  claims,  it  is  shown  in  this  paper  that  the  original  version 
of  the  RNG  K  —  e  model  substantially  underpredicts  the  reattachment  point  in  the  backstop 
problem  -  a  deficiency  that  is  traced  to  the  modeling  of  the  production  of  dissipation  term. 
However,  with  the  most  recent  improvements  in  the  RNG  K  —  e  model  proposed  by  Yakhot 
and  co-workers,  excellent  results  for  the  backstep  problem  are  now  obtained.  Interestingly 
enough,  these  results  are  not  that  sensitive  to  the  details  of  the  near  wall  treatment. 
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1.  INTRODUCTION 


During  the  past  two  decades,  the  K  —  e  model  has  played  a  central  role  in  the  calculation 
of  many  of  the  turbulent  flows  of  scientific  and  engineering  interest.  In  this  two-equation 
model,  the  turbulence  length  and  time  scales  are  built  up  from  the  turbulent  kinetic  energy 
K  and  dissipation  rate  e  which  are  obtained  from  separate  modeled  transport  equations.  The 
Reynolds  stress  tensor  is  represented  by  an  eddy  viscosity  model  constructed  from  length  and 
time  scales  based  on  K  and  e.  It  was  during  the  early  1970’s  that  the  original  version  of  the 
K  —  e  model  was  developed  by  Launder  and  co-workers  [1,  2]  based  on  a  phenomenological 
approach  which  incorporated  empirically  many  of  the  earlier  ideas  of  Kolmogorov  [3].  In  the 
intervening  years,  a  variety  of  modifications  have  been  made  to  the  K  —  e  model  to  account 
for  other  complicating  features  such  as  near-wall  turbulence,  body  forces  and  compressibility 
effects  to  name  a  few  [4-7]. 

One  of  the  major  criticisms  that  has  been  leveled  against  the  K  —  e  model  is  that  it  is 
ad  hoc  -  namely,  it  is  not  derived  from  the  Navier-Stokes  equations  in  a  systematic  fash¬ 
ion.  However,  Yakhot  and  Orszag  [8]  recently  derived  a  version  of  the  K  —  e  model  using 
Renormalization  Group  (RNG)  methods.  In  this  approach,  an  expansion  is  made  about  an 
equilibrium  state  with  known  Gaussian  statistics  by  making  use  of  the  correspondence  prin¬ 
ciple  wher  the  effects  of  mean  strains  are  represented  by  a  random  force.  Bands  of  high 
wavenumbers  (namely,  the  small  scales)  are  systematically  removed  and  space  is  rescaled 
in  a  manner  analogous  to  that  employed  in  the  study  of  phase  transitions.  The  removal  of 
only  the  smallest  scales  gives  rise  to  subgrid  scale  models  for  large-eddy  simulations  whereas 
the  removal  of  successively  larger  scales  ultimately  leads  to  Reynolds  stress  models.  At  high 
turbulence  Reynolds  numbers,  the  RNG  based  K  —  e  model  of  Yakhot  and  Orszag  [8]  is  of 
the  same  general  form  as  the  standard  K  —  e  model;  however,  the  constants  are  calculated 
explicitly  and  assume  somewhat  different  values.  One  major  difference  between  the  RNG 
and  standard  K  —  e  model  lies  in  the  near  wail  treatment.  The  RNG  model  can  be  inte¬ 
grated  directly  to  a  solid  boundary  without  the  need  for  ad  hoc  wall  damping  functions. 
Several  applications  of  this  RNG  K  —  e  model  have  been  subsequently  reported  by  Yakhot, 
Orszag  and  co-workers.  Probably  the  most  notable  one  -  since  it  deals  with  a  complex  shear 
flow  with  separation  -  is  for  turbulent  flow  over  a  backward  facing  step  (see  Karniadakis  et 
al.  [9]).  It  is  well-known  that  the  standard  K  —  e  model  underpredicts  the  reattachment 
point  in  this  problem  -  a  deficiency  that  has  been  widely  discussed  in  the  literature  since 
the  1980/81  Stanford  Conference  on  Complex  Turbulent  Flows  [10].  Karniadakis  r/  al.  [9] 
claimed  that  the  RNG  K  —  e  model  predicts  reattachment  7.3  step  heights  downstream  of 
the  step  corner  for  the  Kim,  Kline  and  Johnston  [11]  test  case  -  a  result  that  is  extremely 
close  to  the  experimental  value.  This  apparent  improvement  in  the  predictions  of  the  RNG 
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based  K  —  e  model  were  attributed  largely  to  the  better  treatment  of  near  wall  turbulcjice 
effects  [9]. 

Subsequent  to  the  publication  of  the  paper  by  Karniadakis  et  al.  [9],  several  adjustments 
have  been  made  in  the  RNG  K  —  e  model.  An  error  in  the  calculation  of  the  constant 
Cci  in  the  dissipation  rate  transport  equation  was  corrected  (see  Smith  and  Reynolds  [12] 
and  Yakhot  and  Smith  [13]).  Furthermore,  a  double  expansion  technique  was  introduced  by 
Yakhot  et  al.  [14]  in  an  effort  to  yield  improved  models  for  the  production  of  dissipation  term 
that  can  better  accommodate  large  strain  rates.  In  the  latter  paper,  the  backstep  problem 
for  the  Kim,  Kline  and  Johnston  [11]  test  case  was  re-computed  yielding  excellent  results  for 
this  latest  version  of  the  RNG  K  —  e  model.  Consequently,  we  feel  that  there  is  a  need  to 
clarify  several  questions: 

(1)  Can  the  original  version  of  the  RNG  K  —  e  model  yield  good  results  for  the  backstep 
problem? 

(2)  Does  the  latest  version  of  the  RNG  K  —  e  model  yield  excellent  results  for  the  backstep 
problem  due  to  the  change  in  the  constant  Cei,  the  addition  of  the  new  model  for  the 
production  of  dissipation  term,  or  some  other  combination  of  factors? 

While  the  first  question  was  answered  to  the  negative  in  a  recent  paper  by  Thangam  and 
Speziale  [15],  there  are  still  some  open  questions  regarding  the  near  wall  modeling  that  need 
to  be  addressed.  In  the  sections  to  follow,  an  effort  will  be  made  to  fully  clarify  these  issues 
and  to  gain  a  better  understanding  of  the  role  of  the  various  modeled  terms  in  the  prediction 
of  turbulent  separated  flows. 


2.  THE  RNG  K  -  e  MODEL  AND  THE  BACKSTEP  PROBLEM 


The  RNG  procedure  of  Yakhot  and  Orszag  [8]  for  incompressible  turbulent  flows  yields 
the  renormalized  equations  of  motion 


du,  _  dui 
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where  Ui  is  the  mean  velocity,  P  is  the  mean  pressure,  u  is  the  molecular  viscosity  and  ut  is 
the  eddy  viscosity.  At  high  turbulence  Reynolds  numbers,  the  eddy  viscosity  takes  the  form 
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are,  respectively,  the  turbulent  kinetic  energy  and  dissipation  rate.  is  a  dimensionless 
constant  which  was  computed  to  be  0.085  using  the  RNG  methodology.  The  turbulent 
kinetic  energy  and  dissipation  rate  are  typically  obtained  from  modeled  transport  equations 
which,  at  high  Reynolds  numbers,  take  the  form 
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where  V  is  the  turbulence  production  given  by 
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is  the  mean  rate  of  strain  tensor.  The  modeled  transport  equations  (5)  -  (6)  contain  four 
coefficients;  Cei,  Ceji  <^k  and  In  the  standard  K  —  e  model,  these  coefficients  are  con¬ 
stants  which  are  obtained  from  benchmark  experiments  for  equilibrium  turbulent  boundary 
layers  and  isotropic  turbulence.  On  the  other  hand,  in  the  original  version  of  the  RNG  K  —  e 
model  of  Yakhot  and  Orszag  [8],  the  coefficients  are  constants  which  are  calculated  explicitly 
by  the  theory.  In  the  newest  version  of  the  RNG  K  —  e  model  developed  by  Yakhot  et  al. 
[14],  corrections  to  the  constants  in  the  dissipation  rate  equation  made  by  Yakhot  and  Smith 
[13]  were  implemented  along  with  a  modification  of  the  production  of  dissipation  term  to 
account  for  large  strain  rates.  These  coefficients  can  be  summarized  as  follows: 


Standard  K  —  e  Model  [2] 

=  0.09,  C,x  =  1.44,  Ce2  =  1.92 
(Tfi  =  1.0,  <7e  =  1.3 

Original  RNG  K  —  e  Model  [8] 

=  0.085,  C,t  =  1.063,  C,2  =  1.72 
(Tk  =  0.7179,  a,  =  0.7179 
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New  RNG  K  -£  Model  [14] 


=  0.085,  Cel  =  1.42  - 


-  V/Vo) 


Ce2  =  1.68,  an  =  0.7179,  =  0.7179 

where  rj  =  SKfe^  S  =  t]q  =  4.38  and  0  =  0.015.  The  new  RNG  K  —  e 

model  is  obtained  by  implementing  the  corrections  of  Yakhot  and  Smith  [13]  (wherein  Cc\ 
was  changed  from  1.063  to  1.42  and  Ce2  from  1.72  to  1.68)  along  with  a  strain-dependent 
modification  to  Cei  which  is  significant  for  flows  with  large  strain  rates.  In  so  far  as  the 
latter  modification  is  concerned,  r}Q  is  the  fixed  point  for  homogeneously  strained  turbulent 
flows  and  is  a  constant  which  was  evaluated  to  yield  a  von  Karman  constant  k  ^  0.4  (see 
Yakhot  et  al.  [14]). 

The  test  case  to  be  considered  is  turbulent  flow  over  a  backward  facing  step  (see  Figure 
1).  Model  predictions  will  be  compared  with  the  experimental  data  of  Kim,  Kline  and 
Johnston  [11]  as  updated  by  Eaton  and  Johnston  [16]  for  the  1980/81  Stanford  Conference 
on  Complex  Turbulent  Flows.  For  this  flow  configuration,  the  expansion  rati"'  E  (step 
height:  outlet  channel  height)  is  1:3  and  the  Reynolds  number  Re  is  132,000  based  on  the 
inlet  centerline  mean  velocity  and  the  outlet  channel  height. 

The  fully-developed  mean  velocity  is  of  the  two-dimensional  form 

u  =  u(x,i/)i-l-t;(a:,y)j.  (12) 


Hence,  the  governing  field  equations  to  be  solved  are  given  by 
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where  the  turbulence  production  V  takes  the  form 


'P  =  2i/x 


(du  dvV 


dv^ 


given  that  the  eddy  viscosity  i/j  is  as  defined  in  (3). 

Equations  (13)  -  (17)  are  solved  subject  to  the  following  boundary  conditions: 

(a)  The  inlet  mean  velocity  profile  u  is  specified  five  step  heights  upstream  of  the  step 
corner.  This  profile  is  matched  to  the  experimental  data  [16]  by  performing  a  separate 
developing  channel  flow  calculation  which  also  provides  companion  profiles  for  K  and 
e  at  the  inlet. 

fb)  Extrapolated  outflow  conditions  are  applied  thirty  step  heights  downstream  of  the  step 
corner. 

(c)  The  law  of  the  wall  is  applied  at  the  upper  and  lower  walls  in  the  three  layer  form  [15]: 

2/+,  for  ?/+  <  5 

u'*'  =  <  —3.05  +  5  y'’",  for  5  <  y'*'  <  30  (19) 

5.5  +  2.5  £n  y"*^,  for  y"^  >  30 

where  y"*"  =  yUrjv,  =  u/ur  and  Ur  is  the  friction  velocity  constructed  from  the  wall 
shear  stress  in  the  standard  manner.  The  law  of  the  wall  (18)  is  used  in  conjunction 
with  the  boundary  condition  of  vanishing  normal  derivative  of  K  at  the  wall  along 
with  the  log  layer  formulas 

K  A'3/2 

4  =  £  =  cj'-—  (20) 

(where  k  %  0.4  is  the  von  Karman  constant)  which  are  suitably  interpolated  to  establish 
boundary  conditions  for  K  and  e  at  the  first  grid  point  away  from  the  wall  [17].  The 
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friction  velocity  Ur  is  approximated  by  its  log  law  form  Ur  in  order  to 

avoid  problems  with  singularities  at  the  reattachment  point  where  the  wall  shear  stress 
vanishes. 


The  equations  of  motion  are  solved  in  their  time  dependent  form  by  a  finite  volume 
method.  A  variable  200  x  100  mesh,  which  concentrates  mesh  points  near  the  step  corner, 
is  used.  This  resolution  is  required  to  ensure  the  numerical  accuracy  of  the  results  (many 
earlier  computations  of  the  backstep  problem  were  under-resolved;  see  Thangam  and  Speziale 

[15]  for  a  discussion  of  this  issue  as  well  as  for  other  details  of  the  numerical  method).  The 
computations  were  conducted  on  a  CRAY-YMP  computer.  Approximately  15  minutes  of 
CPU  time  is  required  to  obtain  a  steady  state  solution  for  a  given  test  case. 

3.  DISCUSSION  OF  THE  RESULTS 

First,  we  will  present  computed  results  for  the  standard  K  —  £  model  since  this  is  an  im¬ 
portant  benchmark  against  which  the  other  results  will  be  judged.  All  model  predictions  will 
be  compared  with  experimental  data  for  the  Kim,  Kline  and  Johnston  test  case  as  updated  by 
Eaton  and  Johnston  [16]  for  the  1980/81  Stanford  Conference  on  Complex  Turbulent  Flows. 
In  Figures  2(a)-(b)  the  mean  velocity  streamlines  and  mean  velocity  profiles  obtained  from 
the  standard  K  —  e  model  are  compared  with  the  experimental  data.  The  standard  K  —  e 
model  predicts  reattachment  at  Xr! H  «  6.25  -  a  value  that  constitutes  approximately  a 
12%  underprediction  of  the  mean  experimental  reattachment  point  of  Xr! H  w  7.1.  Earlier 
reported  results  [10]  of  a  20-25%  underprediction  of  the  reattachment  point  appear  now  to 
have  been  in  error  due  to  insufficient  numerical  resolution  as  discussed  by  Thangam  and 
Speziale  [15]  and  Avva,  Kline  and  Ferziger  [18]. 

One  of  the  criticisms  that  h^ls  been  leveled  against  most  of  the  K  —  e  model  calculations 
for  the  backstep  problem  lies  in  the  use  of  wall  functions.  It  is  well-known  that  wall  functions 
do  not  formally  apply  to  separated  flows.  However,  since  the  separation  point  is  fixed  -  and, 
perhaps,  since  a  large  proportion  of  the  turbulent  kinetic  energy  is  not  associated  with  the 
separated  zone  -  it  now  appears  that  wall  functions  do  not  give  rise  to  major  errors  in 
the  backstep  problem,  unlike  in  other  separated  flows.  A  systematic  study  of  near  wall 
turbulence  models  by  the  authors  has  indicated  that  they  generally  do  not  give  rise  to  more 
than  a  5%  improvement  in  the  results  for  the  backstep  problem  obtained  from  well-resolved 
computations  using  wall  functions.  In  Table  1,  the  reattachment  point  prediction  of  the 
standard  K  —  e  model  using  three-layer  wall  functions  is  compared  with  experimental  data 

[16]  as  well  as  with  the  result  obtained  using  the  £isymptotically  consistent  near  wall  model 
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of  Speziale,  Abid  and  Anderson  [19].  Here,  this  particular  near  wall  model  is  chosen  since  it 
was  shown  by  those  authors  to  be  better  behaved  than  four  independent  near-wall  turbulence 
models  that  were  recently  reviewed  by  Patel,  Rodi  and  Scheuerer  [20]  for  boundary  layer 
flows.  It  is  clear  from  Table  1  that  the  use  of  an  asymptotically  consistent  near  wall  model 
does  not  lead  to  a  significant  improvement  in  the  results  obtained  using  wall  functions. 
Furthermore,  it  should  be  noted  that  the  use  of  an  asymptotically  inconsistent  near  wall 
model  (some  of  which  are  discussed  in  Patel,  Rodi  and  Scheuerer  [20])  can  lead  to  even 
worse  results  than  those  obtained  using  wall  functions!  Consequently,  we  believe  that  it  is 
acceptable  to  use  wall  functions  when  comparing  the  predictions  of  a  variety  of  turbulence 
models  in  the  backstep  problem. 

Now  we  will  present  computed  results  for  the  original  version  of  the  RNG  K  —  e  model 
(Yakhot  and  Orszag  [8]).  In  Figures  3(a)  -  (b)  the  computed  mean  velocity  streamlines 
and  mean  velocity  profiles  are  compared  with  the  experimental  data  of  Eaton  and  Johnston 
[16].  The  surprising  finding  is  that  the  predicted  mean  reattachment  point  of  X  h/H  «  4.0 
constitutes  over  a  40%  underprediction  of  the  mean  experimental  result  of  XnfH  «  7.1. 
This  deficiency  -  as  well  as  a  plausible  explanation  for  its  cause  -  was  pointed  out  recently 
by  Thangam  and  Speziale  [15].  In  a  homogeneous  shear  flow,  the  eddy  viscosity  grows 
exponentially: 

vt  ~  exp(A<*) 


where  t”  is  the  time  (nondimensionalized  by  the  shear  rate)  and  A  is  the  growth  rate  given 
by  (see  Speziale  [21]) 


A  = 


1/2 


(21) 


which  becomes  singular  for  C^i  =  1.  Consequently,  for  Cji  =  1.063,  the  growth  rate  of  the 
eddy  viscosity  will  be  overly  large.  When  a  model  overpredicts  the  eddy  viscosity  it  will 
be  too  dissipative  -  a  deficiency  that  will  lead  to  an  underprediction  of  the  separated  flow 
region  in  the  backstep  problem.  Earlier  reported  results  by  Karniadakis  et  al.  [9]  for  the 
backstep  problem  -  which  suggested  that  the  model  predictions  were  excellent  -  now  appear 
to  have  been  in  error. 

With  the  correction  of  C\i  to  1.42  by  Yakhot  and  Smith  [13],  the  problem  of  overpredicting 
the  growth  rate  of  homogeneous  shear  flow  is  eliminated.  Speziale,  Gatski  and  Fitzmaurice 
[22]  showed  that  with  the  new  value  of  C^i,  the  RNG  K  —  e  model  yields  excellent  results 
for  homogeneous  shear  flow.  However,  an  additional  production  term  in  the  dissipation 
rate  transport  equation  was  uncovered  which  Yakhot  and  Smith  [13]  were  not  able  to  close 
systematically.  The  new  term  was  shown  to  be  important  for  large  strain  rates.  If  this 
term  is  neglected,  a  revised  RNG  K  —  e  model  is  obtained  which  is  of  the  same  general 
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form  as  the  Yakhot  and  Orszag  [8]  model;  only  the  values  of  the  constants  are  altered.  In 
Figures  4(a)-(b),  the  mean  velocity  streamlines  and  mean  velocity  profiles  obtained  from  the 
Yakhot-Smith  RNG  K  —  e  model  (ignoring  the  additional  production  term)  are  shown.  It 
is  clear  that  the  model  drastically  overpredicts  the  reattachment  point  (i.e.,  XrJ H  ?=:  9.7  in 
comparison  to  the  mean  experimental  result  of  7.1).  The  origin  of  this  erroneous  prediction 
is  clear:  with  the  additional  production  term  in  the  e-transport  equation  neglected,  the 
Yakhot-Smith  model  yields  a  von  Karman  constant  of  0.23  instead  of  the  traditional  value 
of  0.4  (see  Yakhot  and  Smith  [13]).  This  strongly  indicates  that  the  additional  production 
term  cannot  be  neglected  in  wall  bounded  turbulent  flows. 

Yakhot  et  al.  [14]  recently  developed  a  model  for  the  additional  production  term  in  the 
£-transport  equation  by  means  of  a  scale  expansion.  More  precisely,  an  expansion  of  this 
term  was  made  in  the  ratio  of  the  turbulent  to  mean  strain  time  scale  //  =  SK/e,  where 
S  =  (25,j5ij)^^^  is  the  norm  of  the  mean  rate  of  strain  tensor.  The  interesting  finding  was 
that  no  finite  truncation  of  the  expansion  for  the  production  of  dissipation  term  in  powers  of 
1]  suffices  -  terms  of  all  orders  must  be  retained  to  satisfy  the  crucial  weak  and  strong  strain 
limits  [14].  This  complication  eliminated  the  possibility  of  obtaining  a  closed  form  solution 
for  this  term.  A  Fade  approximation  was  made  which  had  one  undetermined  constant  /?  that 
was  evaluated  empirically  by  setting  the  von  Karman  constant  to  0.4.  This  leads  to  a  net 
production  of  dissipation  term  of  the  traditional  form  Cii{e/K)V  with  a  variable  coefficient 
given  by 

^  ...  V{l-ri/Vo) 

where  f3  =  0.015  and  t/o  =  4.38  is  the  fixed  point  for  equilibrium  homogeneous  flows.  The 
addition  of  the  second  term  to  the  r.h.s.  of  (22)  -  which  becomes  signihcant  for  large  strain 
rates  J]  >  Tjo  -  is  the  feature  that  distinguishes  the  Yakhot  et  al.  [14]  model  from  the  Yakhot 
and  Smith  [13]  model. 

In  Figures  5(a)-(b),  the  computed  mean  velocity  streamlines  and  mean  velocity  profdes 
for  the  Yakhot  et  al.  [14]  RNG  K  —  e  model  are  compared  with  experimental  data  [16].  The 
model  predicts  reattachment  at  Xr/H  ^  6.7;  a  result  that  is  within  5%  of  the  experimental 
mean  reattachment  point  of  7.1.  Most  notably,  this  new  RNG  K  —  e  model  -  unlike  all  of 
the  other  models  -  yields  a  significant  secondary  separation  bubble  below  the  corner  of  the 
step  that  is  more  in  line  with  the  experimental  data  [23].  Yakhot  et  al.  [14]  showed  that 
when  this  model  is  extended  to  include  an  anisotropic  eddy  viscosity,  the  predicted  mean 
reattachment  point  is  almost  identical  to  the  experimental  value  of  XrJ H  ^  7.1. 


8 


4.  CONCLUDING  REMARKS 


Three  independent  versions  of  the  RNG  K  —  e  model  have  been  applied  to  the  test  case 
of  turbulent  flow  over  a  backward  facing  step  in  order  to  assess  its  performance  in  separated 
flows.  The  main  objective  of  this  study  was  to  better  understand  the  impact  that  the  recent 
changes  in  the  Yakhot  h  Orszag  RNG  K  —  e  model  have  had  on  the  predictive  capabilities  of 
the  model  -  particularly  in  a  complex  shear  flow  of  engineering  interest.  It  was  found  that  the 
original  version  of  the  RNG  K  —  e  model  substantially  underpredicts  the  reattachment  point 
for  the  Kim,  Kline  and  Johnston  test  case  due  to  an  overprediction  of  the  eddy  viscosity 
(the  model  is  far  too  dissipative).  However,  the  latest  version  of  the  RNG  K  —  c  model 
[14]  -  which  is  an  extension  of  the  Yakhot  and  Smith  [13]  model  to  accommodate  large 
strain  rates  in  the  production  of  dissipation  term  -  yields  excellent  results  for  the  backstep 
problem.  The  mean  reattachment  point  is  predicted  to  within  5%  of  the  experimental  value 
in  comparison  to  the  standard  K  —  e  model  which  yields  a  12%  error.  Furthermore,  a 
secondary  separation  bubble  was  obtained  that  is  much  more  in  agreement  with  that  which 
is  observed  in  experiments.  These  results  are  strongly  dependent  on  the  presence  of  the 
variable  strain  rate  term  in  (21)  and  the  calculations  vividly  demonstrate  the  sensitivity  of 
the  reattachment  point  to  the  values  of  Cf\  and  (7^2-  The  encouraging  point,  however,  is  that 
the  results  are  not  that  sensitive  to  the  details  of  the  near  wall  treatment  and  are  obtained 
with  no  ad  hoc  adjustments  of  the  constants.  Based  on  these  results,  we  feel  that  the  current 
version  of  the  RNG  K  —  £  model  can  be  a  useful  turbulence  model  for  practical  engineering 
and  scientific  calculations. 
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NEAR-WALL  MODEL 

XrIH 

Three-layer  wall  functions 

6.25 

Speziale  et  al.  [19] 

6.40 

Experimental  data  [16] 

7.1 

Table  1.  The  reattachment  point  in  the  backstep  problem:  Comparison  of  the  predictions  of 
different  near  wall  corrections  to  the  standard  K  —  e  model  with  experimental  data. 
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Figure  1.  Schematic  of  turbulent  flow  over  a  backward  facing  step. 
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(b)  Dimensionless  Mean  Velocity  Profiles 


Figure  2.  Comparison  of  the  predictions  of  the  standard  K  —  5  model  [2]  with  experimental 
data  [16].  (a)  Streamlines  and  (b)  mean  velocity  profiles. 
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Figure  3.  Comparison  of  the  predictions  of  the  RNG  K  —  e  model  of  Yakhot  and  Orszag  [8] 
with  experimental  data  [16].  (a)  Strezunlines  and  (b)  mean  velocity  profiles. 
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(b)  Dimensionless  Mean  Velocity  Profiles 

Figure  4.  Comparison  of  the  predictions  of  the  RNG  K  —  £  model  of  Yakhot  and  Smith  [13] 
-  where  the  additional  production  of  dissipation  term  is  neglected  -  with  experimental  data 
[16].  (a)  Streamlines  and  (b)  mean  velocity  profiles. 
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Figure  5.  Comparison  of  the  predictions  of  the  RNG  K  -e  model  of  Yakhot  et  al.  [14]  with 
experimental  data  [16].  (a)  Streamlines  and  (b)  mean  velocity  profiles. 
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